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ADRENERGIC RECEPTOR SNP FOR IMPROVED MILKING CHARACTERISTICS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] Not applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 

[0002] NotappUcable. 

BACKGROUND OF THE INVENTION 
[0003] The present invention relates, in general, to the field of molecular biology, human 

and bovine genetics, and desirable milking characteristics. In particular this invention provides 
a method for using genetics to predict the milking characteristics of cows. 
[0004] Various pubUcations or patent are referred to throughout this application to 

describe the state of the art to which the invention pertains. Each of this publications or patents 
is incorporated by reference herein. Complete citations of scientific publications are set forth in 
the text or at the end of the specification. 

[0005] Domestic animals are bred for improvCTient to the usefiil attributes of the 

animals. For cows, one very usefijl attribute is their milking characteristics, both volume of 
milk and rate of milking. Many dairy cows are artificially inseminated using semen fipom sires 
purchased firom companies engaged in selling quality improvement genetics for milking cows. 
However, while the overall statistical average of the milking characteristics of the daughters of a 
particular sire can be measured, it is yet not possible to predict Avith any degree of success the 
milking characteristics of a particular progeny created by breeding a cow with a particular bull. 
Better understanding of the genetics of milking traits will be needed in order to more accurately 
predict the milking characteristics of a particular daughter. 

[0006] It is known that some milking characteristics are linked to other traits. A 

relationships between adrenoceptor concentrations and milkability in primiparous cows have 
been elucidated by Roets, et al.(^elationship between milkability and adrenocq)tor 
concentrations in teat tissue in primiparous cows" J Dairy Sci 69(12):2131-3131, 1986). Later, 
Roets et al. ('^Relationship between numbers of alpha2- and beta2-adrenoceptors in teat tissue 
and blood cells and milkability of primiparous cows", J Dairy Sci 72(12):3304-33-13, 1989) 
reported on flie numbers of alpha2' and beta2-adrOToceptors in teat tissue and blood cells and 

-1- 



wo 2005/030789 



PCTAJS2004/030774 



milkability in primiparous cows. Subsequently, Roets et al.CTRelationship between numbers of 
alpha2- and beta2-adrenoceptors on blood cells of bulls and milkability of their daughters", J 
Dairy Sci, 62 567-575, 1995) analyzed the numbers of alpha2- and beta2-adrenoceptors on the 
blood cells of bulls and the milkability of their daughters, hi addition, Blum et 
al.("Catecholamines, oxytocen and milk removal in dairy cows**, J Dairy Res, 56(2): 167-1 77, 
1989) reported tibat treatment with phentolamine (an alpha-adrraergic blocker) and isoproterenol 
(a beta-adrmoceptor agonist) each faciUtated milk removal. Brown et al. ('^Relationship of 
milking rate to somatic cell count", J Dairy Sci, 69(3):850-854, 1986) raised the possibility that 
milking speed was correlated with botihi somatic cell score (SCS) and mastitis. Many studies 
have established correlations between mastitis and SCS (see, for example, MacMillan et al., 
"Associations between dry cow therapy, clinical mastitis, and somatic cell count score with milk 
and fat production in ten New Zealand dairy herds*', J Dairy Sci, 66(2):259-265, 1983; Reneau, 
"Effective use of dairy herd improvement somatic cell counts in mastitis control", J Dairy Sci, 
69(6): 1708-1 720, 1986; and McClelland, "A comparison of objective and subject measures of 
milking speed on Canadian Holstein-Friesians", University of Guelph, 1983). Further, a 
positive correlation between milking speed and SCS was confirmed by Boettcher et al. 
('Development of an udder health index for sire selection based on somatic cell score, udder 
conformation, and milkmg speed", J Dairy Sci, 81(4):1157-1168, 1998). 

[0007] The DNA sequence in and around the bovine beta2-adrenergic recqptor (ADRB2) 

gene was in public genome databases (accession numbers Z86037; Einspanier et al., 
"Expression of the beta2 adrenergic receptor in the cattle", direct submission to Genbank, 
accession No. Z86037, 1997 and AF331034; Schimpf et al., "Genetic mapping of the ADRB2 
gene on cattle chromosome T\ Anim. Genet. 32(6):390, 2001). These sequences included the 
ADRB2 coding region (1257 bases inclusive from the ATG start to the TAA stop codon), as 
well as the 223 bases upstream from the ATG start codon and the 550 bases downstream from 
the TAA stop codon. 

[0008] What is needed is a way of determining which bidls will produce cows with 

improved milkability and a method of improving the results of a breeding program for 
milkability. 



BRIEF SUMMARY OF THE INVENTION 
[0009] Disclosed herein is a method for breeding cows for improved milking 

characteristics including the step of screening the genotype of the parents of the cow for the 
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allele associated with a desired SCS phenotype. The method has the steps of obtaining a DNA 
sample from a bull to be tested for the desired SCS phenotype; and detecting the presence of an 
adenine at position 11 in a gene encoding a bovine beta2-adrenorecepton The method can be 
performed by direct sequencing, primer extension, restriction length fragment polymorphism, 
and allele-specific hybridization. 

[00010] In another embodiment, a method of identifies a bull whose daughter cows have a 
short milking duration This method has fte steps of obtaining a sample of DNA from a bull; 
combining the DNA with a pair of PGR probes comprising SEQ IDs 1 and 2 or SEQ IDS 3 and 
4; incubating the DNA under conditions pemiitting the DNA bounded by the PGR probes to 
produce DNA ampKcons; isolating the DNA amplicons; combining the DNA amplicons with a 
restriction enzyme specific for CCCGGG for a sufficient time to produce a mixture of DNA 
fragments from the amplicons comprising CCGGGG; applying the DNA fragment mixture to a 
gel and permitting migration of the mixture components for a time sufficient for them to 
separate; and observing the sizes of DNA on the gel, with the largest fragments being correlated 
with the A genotype and with better SCS phenotype, and the smaller fragments being associated 
with the C genotype and less desirable SCS phenotype. 

[00011] In another embodiment a milking attribute PCR-RFLP kit contains a pair of . 
primers which flank the 1 1*^ nucleotide position of the bovine beta2"adrenoreceptor gene, and a 
restriction enzyme specific for the CCCGGG site. The restriction enzyme can be Smal. The 
primer pairs are selected fit>m pair 1 (SEQ ID Nos 1 and 2) or from pair 2 (SEQ ID Nos 3 and 
4). 

[00012] Other objects, features and advantage of the present invention will become 
apparent to one of skill in the art after review of the specification and claims. 

BRIEF DESCRIPTION OF TEDE SEVERAL VIEWS OF THE DRAAVINGS 
[00013] Not applicable. 

DETAILED DESCRIPTION OF THE DSTVENTION 
[00014] It is disclosed here that genetic differences in the allele for the beta2-adrenoceptor 
(ADR2) gene are associated with the milking characteristics of milking cows. This makes 
possible the use of genetic tools and analysis in breeding of cows to start to predict with 
accuracy tiie milking characteristics of daughters of specific bulls bred to specific maternal 
cows. 
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[0001 5] We began the work described here by proposing that DNA sequence 
polymorphisms in and aronnd bovine beta2-adrenoceptor (ADRB2) may cause or be associated 
with genes that cause the observed associations. First, we subcloned and sequenced bovine 
ADRB2, extending the region of knovm sequence beyond what was previously publicly 
available. Sequence obtained through this work, not previously in the public domain, includes 
876 additional bases upstream from that described above. These bases are a total of 1099 that 
are S' to the ATG start codon. An additional 695 bases downstream fiom the previously 
reported bases also were sequenced. These bases include a total of 1245 bases 3' to the TAA 
stop codon. 

[0001 61 With this information in hand, we were then able to design oligonucleotide 
primers to be used to amplify and perform locus-specific re-sequencing of genomic DNA 
spanning regions of the bovine ADRB2 gene. These primers were used to amplify genomic 
DNA from 24 Holstein and 12 Brown Swiss dairy cows. After comparative genomic sequence 
was obtained, we analyzed the sequence using polyPhred software (Washington University, St. 
Louis, MO) and identified, among others, a Single Nucleotide Polymorphism (SNP) sequence at 
position 1 1 , inclusive of the start ATG (hereafter referred to as "Al IC). The published 
sequOTce indicates the presence of a cytosine (C) nucleotide. We discovered a variant allele 
with an adenine (A) nucleotide at that position. Substitution of A for C at the nucleotide 
position 1 1 altered the amino acid sequence of ADRB2 protein by replacing a proline with a 
histidine at amino acid position 4 of the protein. 

[000171 Subsequently we designed and validated a polymerase chain reaction-restriction 
fragment length polymorphism (PCR-RPLP) assay able to detect the presence of A and/or C 
nucleotide(s) at the Al IC locus. This assay is based on the presence or absence of the smal 
restriction enzyme recognition site (CCCGGG) in DNA fragments which have been amplified 
by one of two primer pairs. The A allele at the Al 1 C locus changes this sequence to CACGGG, 
which is not recognized, therefore is not digested hy smal and consequently is identified as a 
heavier nucleotide in the PCR-RFLP test. 

[00018] This polymorphism was found to be associated with miUdng characteristics. In 
particular, the A allele at the AUG locus is associated with higher somatic cell score (SMS) and 
therefore with a faster milking speed. This makes is possible to breed dairy cattle to increase the 
proportion of the herds having the A allele at the Al IC locus, to improve flie overall milking 
speed of the herd. The A allele is effective in both the heterozygous and homozygous conditions 
to improve the performance of the animal. 
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[00019] 'Tolymorphism" refers generally to the ability of an organism or gene to occur in 
two or more different forms. In particular for purposes of the present invention, 
"polymorphism" refers to two or more different fomis of the same gene. 
[00020] "Single Nucleotide Polymorphism" or "SNP" ref^s to a polymorphism that 
results from a difference in a single nucleotide. 

[00021] A "Restriction Enzyme" refers to an endonuclease which bins to double stranded 
DNA at a specijBc nucleotide sequence and then, if both strands of the DNA lack the appropriate 
modification at that sequence (including but not limited to methylation), cleaves the DNA either 
at the recognition sequence or at anoth^ site in the DNA molecule. Restriction razymes are 
denoted by three-letter abbreviations followed by a strain designation and/or a Roman numeral 
distinguishing different enzymes from the same species or strain. Recognition sequences are 
written 5* to 3' for one strand only. Examples of restriction enzymes include Smal, BamHi, 
Bell, EcoRI, Hindm, and Xbal. 

[00022] The term "Allele" refers generally to any of one or more alternative forms of a 
given gene or DNA segment; both or all alleles of a given gene are concerned with the same trait 
or characteristic, but the product or function coded for by a particular allele differs qualitatively 
and/or quantitatively from that coded for by other alleles of that gene. In a diploid cell or 
organism the members of an allelic pair (i.e., the two alleles of a given gene) occupy 
corresponding positions (loci) on a pair of homologous chromosomes; if these alleles are 
genetically identical the cell or organism is said to be homozygous. If flie alleles are genetically 
different, the cell or organism is said to be heterozygous with respect to flie particular gene. 
[00023] The polynucleotides of the present invention may be prepared by two general 
methods: (1) They may be synthesized from appropriate nucleotide triphosphates, or (2) they 
may be isolated from biological sources. Both methods utilize protocols well known in the art. 
The availability of nucleotide sequence information enables preparation of an isolated nucleic 
acid molecule of the invention by oligonucleotide synthesis. Synthetic oligonucleotides may be 
prepared by the phosphoramidite method employed in the Applied Biosystems 37A DNA 
Synthesizer or similar devices. The resultant constmct may be purified according to methods 
known in the art, such as high performance lic^uid chromatography (HPLC). Complementary 
segments thus produced may be annealed such that each segment possesses appropriate cohesive 
termini for attachment of an adjacent segment. Adjacent segments maybe ligated by annealing 
cohesive termini in the presence of DNA ligase to construct an entire long double-stranded 
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molecule. A synthetic DNA molecule so constructed may then be cloned and amplified in an 
appropriate vector. 

[00024] The invention will be bett^ understood upon consideration of the following non- 
limiting Examples. Unless ottierwise specified, general cloning procedures, such as those set 
forth in Sambrook et al.. Molecular Cloning, Cold Spring Harbor Laboratory (1989) or 
Ausubel et al. CURRENT Protocols in Molecular Biology, John Wiley & Sons (2001), are 
used. 

Example 1 . Assessment of mKNA levels in leukocytes of cows with fast vs. slow-milking rates. 
[00025] Eight animals were selected by the duration of their milking and grouped as 
either slow or fast. For initial determinations, slow-milking animals were numbered 90, 264, 
279 and 321. Fast-milking animals were numbered 272, 273, 281 and 289, Their blood (about 
15 mL) was collected in EDTA-fiUed Vacutainers (BD, Franklin Lakes, NJ). Known sequences 
jfrom NCBI along with the program "Oligo" were used to design primer sets specific for 
GAPDH, alpha2-, beta 1-, beta2-adrenoreceptors, and beta-arrestin. 

[00026] From the blood samples, RNA was isolated using the RNeasy RNA isolation 
protocol (Qiagen, Valencia, CA). One volume of whole blood was mixed with 5 volumes of 
erythrocyte lysis buffer. The mixture was incubated on ice, and during the incubation mixed by 
vortexing briefly twice. The mixture was next centrifuged at 400 g for 10 min at 4**C. The 
supernatant was discarded and the leukocyte pellet saved. The leukocyte pellet was added to 
more Buffer EL and vortexed to resuspend the leukocytes. The mixture was again centrifuged, 
after which the supernatant was discarded. The leukocyte pellet was loosened by flicking the 
tube. Buffer RLT (Qiagen)(with 10 ^1 jS-mercaptoethanol added per 1 ml buffer) was added to 
the leukocytes; for a starting blood volume of 4.0 ml or less, 2.0 ml was added; and for a greater 
blood volume, 4.0 ml was added. Next the mixture was homogenized until the sample was 
homogeneous. An equal volume of 70% alcohol was added to ttie homogenized lysate and 
vigorously shaken. Finally, the sample, including any precipitate that may have fomied, was 
applied to an RNeasy midi spin column placed in a IS-ml centrifuge tube and centrifuged for 
3000 to 5000 g for 5 min. For samples in excess of 4.0 ml, the extra sample was applied 
successively to tibie same colunm and the flow-dirough discarded. After the flow-through had 
been discarded, 4.0 ml of Buffer RWl was added to tihe RNeasy colunm and centrifuged at 
3000-5000 for 5 min to wash the column. The flow-through was discarded. Next 2.5 ml of 
Buffer RPE (diluted with 4 volumes of ethanol (96-100%) was added to the colunm which was 
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again centrifiiged at 3000-5000 for 2 min. Another 2.5 ml Buffer RPE was added to the column 
which was centrifiiged at 3000-5000 for 5 min to dry the spin-colimm memhrane. KNA was 
eluted from the column adding 150 \xl or 250 \x\ RNase-free water to the column, which was left 
to stand for 1 min and then centrifiiged at 3000-5000 g for 3 min. A second volume of RNase- 
free water was added to the column and spun down. The KNA was then salt-precipitated 
according to the manui^cturer's protocol. 

[00027] Next RNA was treated with DNA-free DNase to remove trace levels of DNA, 
according to flie manufacturer's protocol (Ambion, Austin, TX). To the KNA sample purified as 
described above, were added 0.1 volume of lOX DNase I Buffer and 1 jxl of DNase I (2 units) to 
the RNA, which solution was mixed gently and incubated at 37°C for 20-30 min. Ambion 
DNase Inactivation reagent was resuspended by flicking or vortexing. From that tube, the 
greater of 0.1 volume or 5 ^il was added to the sample and mixed well. The tube was incubated 
for 2 min at room temperature. The tube was next centrifiiged at 10,000 g for about 1 min to 
pellet the DNase inactivation reagent. 

[00028] Spectrophotometry was then used to determine KNA concentration, and 1 fig 
RNA was loaded on a 1% agarose gel using ethidium bromide to ensure quality. 
[00029] Following the manufacturer's protocol, cDNA templates were made from the 
RNA using Qmniscript (Qiagen) for Real Time PGR analysis. The Qmniscript lOX Buffer RT, 
dNTP mix and RNase-free water were first thawed and mixed by vortexing. Qiagen RNase 
inhibitor was diluted to a final concentration of 10 imits/|Lil in Ix Buffer RT and mixed by 
vortexing briefly. Master mix (2.0 ^il of 1 Ox Buffer RT, 2.0 jil dNTP mix, 2.0 jil oligo dT 
primer (10 uM), 0.5 \x\ Rnase inhibitor (10 units/jil), 1.0 jil Omniscript reverse transcriptase and 
double distilled water) was prepared on ice. Master solution was distributed to the various 
tubes, and template RNA (0.6 to 0.68 |il/tube) was added and mixed. The resulting solutions 
were incubated at 37^C for 60 min. 

[00030] All samples were run by LightCycler (Idaho Technology, Inc., Salt Lake City, 
UT) protocol using the SYBR-green kit from Roche Molecular Biochemicals, except for beta2- 
adrenoreceptor, for which the SYBR-green kit (Qiagen) was used. Samples were run in 
diqjlicate, and all were normalized using GAPDH. Cycle thresholds (determined by the second 
derivative growth curve) of samples were then compared between the groups of slow and fast 
milkers. To determine if there were any significant differences between slow and fast milkers, 
the SAS statistical package was used to perform a Pearson correlation coefficient. The results 
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suggested that duration of milking and rate of milk removal are associated with the amount of 
beta2-adrenoreceptor mKNA levels in blood, and it was decided to more thoroughly characterize 
the beta2-adrenoreceptor gene in cattle. 

[00031] It was discovered that there was a polymoiphism of the beta2-adrenoreceptor at 
the 1 1^ nucleotide (A replaced C in some instances, resulting in a substitution of histidine for 
proline at the fourth amino acid). 

Example 2. Associations Between Al IC Genotype and SCS and/or Milking Speed in Dairy 
Cows 

[00032] Bovine DNAs are available from the Cooperative Dairy DNA Repository 
(CDDR) population (Gene Evaluation and Mapping Laboratory, Bldg. 200 Rm 2A, ARS- 
USDA, BARC-East Beltsville, MD 20705). SCS phenotypes have been obtained for all CDDR 
animals. For a subset of the CDDR animals, there also are data on the milking speed (MS). 
First, CDDR animals with MS data were genotyped, along with a subset of the remainder of the 
CDDR animals representing "high" and "low" Somatic Cell Score (SCS) phenotypic classes, 
with the All C assay. DNA samples were obtained and assayed. Data sets were distributed and 
analyzed in duplicate. The analysis revealed a correlation between Al IC and SCS. 
[00033] Six hundred sixty three animals from the CDDR wore genotyped for the Al 1 C 
locus using the either of the following pairs of primers: 

TGGAACTGGCTGAACTGACA (SEQ ID NO 1) 

AGTTGATGGCTTCCTTGTGG (SEQ ID NO 2) 

AGGTCCGCTCGCTGAGG (SEQ ID NO 3) 

GTTCCAGCGTGACGTTTTG (SEQ ID NO 4) 
[00034] Table 1 shows a disproportionate distribution of Al 1 C genotypes into "High" and 
"Low" SCS phenotypic categories. 

Table 1. AllC Genotypes and SCS distribution of 663 CDDR bulls. 





AA 


AC 


CC 


Total 


High # 


9 


88 


248 


345 


Low # 


5 


55 


258 


318 


Total 


14 


143 


506 


663 


Prop . High 


0.643 


0.615 


0.490 


0.520 


Prop . Low 


0.357 


0.385 


0.510 


0.480 
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[00035] In addition. Table 2 contains average estimated transmitting abilities (ETA) for 
SCS and frequencies by Al IC genotype. Results from statistical analyses of the data indicate 
tbat genotypes AA, AC and CC are significantly dififerent from the mean SCS (p<0.025) and 
that genotypes AA and AC are significantly higher than the mean SCS (p<0.01). Furthennore, 
animals with the CC genotype have significantly different (i.e., lower) SCS phenotypes than 
animals with AA and AC genotypes (p<0.05). 

Table 2. Average SCS ETA and frequency by Al IC graotype among 663 CDDR bulls. 

AA AC CC Total 

SCS ETA 314.5 312.8 307.0 0.520 

Genotype Freq. 0.021 0.216 0.763 

[00036] Table 3 shows the distribution of milking speed (MS) classifications according to 
genotype, and Table 4 illustrates the average ETA fr>r MS by genotype. Because of the small 
number of AA animals, there is not a significant association between MS and Al IC genotype at 
the 5% level. However, the trend in Table 3 clearly shows that the CC genotype has a lower 
proportion of bulls transmitting fast MS. 

Table 3. Al IC Genotypes and Milking Speed distribution of 663 CDDR bulls. 





AA 




AC 


CC 




Total 


Fast # 


13 




46 


108 




167 


Slow # 


7 




36 


120 




163 


Total 


20 




82 


228 




330 


Prop. Fast 


0.650 


0 


.561 


0.474 


0 


.506 


Prop . Slow 


0.350 


0 


.439 


0.526 


0 


.494 



[00037] In Table 4, the AA genotype tends to have a higher av^age ETA for MS, 
Table 4. Average ETA for milking speed and frequency by Al IC genotype among 663 CDDR 
bulls. 

AA AC CC 

MS ETA 70 68.81707 64.02609 

Genotype Freq. 0.061 0.248 0.691 

[00038] From these tables it appears that the AA and AC genotypes (that are significantly 
associated SCS) are also associated with higher milking speed. This is in agreanent with 
studies (McClelland, "A comparison of objective and subject m^isures of milking speed on 
Canadian Holstein-Friesians", University of Guelph, 1983, and Boettcher et al., **Development 
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of an udder health index for sire selection based on somatic cell score, udder conformation, and 
milking speed", J Dairy Sci 8 1(4): 11 57-1 168, 1998) that have shown that higher milking speed 
is genetically correlated with higher SCS and mastitis incidence. Since the bulls that had MS 
ETAs represented only a subset (330) of the total group (663), it is likely that with a larger 
sample size, significant association could be shown between MS and tiie Al IC locus. 

[000391 This discovery enables for the first time the testing and selection of animals for 
breeding and/or alternative management practices based on the results of the inventive assay. 
The assay results correlate with the SCS and enable the choice of animals with improved MS 
and mastitis resistance. In addition to using these discoveries to select animals for phenotypes 
related to SCS, these discoveries can also be used to select animals in an attempt to effect 
population changes in MS and mastitis resistance. Applying the assay will enable the selection 
of animals with SCS which will ultimately improve the herd averages for MS and mastitis 
resistance. 

[00040] This invention also includes a kit containing reagents that can be used to identify 
allelic composition at the loci described herein. 

[00041] The present invention is not limited to the embodiments described and 
exemplified above, but is capable of variation and modification within the scope of the appended 
claims. 
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